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Abstract 

A tile Hamiltonian (TH) replaces the actual atomic interactions in a qua- 
sicrystal with effective interactions between and within tiles. We studied 
Al-Co-Cu decagonal quasicrystals described as decorated Hexagon-Boat-Star 
(HBS) tiles using ab-initio methods. The dominant term in the TH counts 
the number of H, B and S tiles. Phason flips that replace an HS pair with a 
BB pair lower the energy. In Penrose tilings, quasiperiodicity is forced by ar- 
row matching rules on tile edges. The edge arrow orientation in our model of 
AICoCu is due to Co/Cu chemical ordering. Tile edges meet in vertices with 
72° or 144° angles. We find strong interactions between edge orientations at 
72° vertices that force a type of matching rule. Interactions at 144° vertices 
are somewhat weaker. 
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I. INTRODUCTION 



Both quasicrystals and ordinary crystals are made of elementary building blocks. In 
crystals, copies of a single building block (known as a unit cell) are arranged side by side to 
cover the space periodically. In quasicrystals, building blocks are arranged to cover the space 
quasiperiodically. Two approaches to build quasilattices have been proposed. One approach 
uses a single unit cell, but allows adjacent cells to overlap. Gummelt M proved that using a 
limited number of overlapping positions between decagons produces a quasicrystal structure. 
Steinhardt and Jeong further proved that the overlapping conditions can be relaxed when 
supplemented by the maximization of a specific cluster density to produce quasilattices. This 
is the "quasi-unit cell" approach. 

In the other "tiling" approach, space is covered with building blocks called tiles. Two or 
more tile types are used 0. No overlapping is allowed, and depending on the way the tiles 
are arranged, quasicrystal structures can be produced. Matching rules between tiles govern 
the local tile configurations by allowing only a subset of all possible arrangements. Globally, 
the matching rules enforce quasiperiodicity Penrose proposed his famous matching 

rules before the discovery of quasicrystalline materials. In 2D, Penrose tiles are fat and 
thin rhombi (Fig |I]a). Edges are assigned arrow decorations which must match for common 
edges in adjacent tiles. In perfect quasicrystals these rules are obeyed everywhere. The very 
restrictive Penrose matching rules are sufficient, but not necessary, to force quasiperiodicity. 
Socolar |J showed that weaker matching rules can still force quasilattices. The less restric- 
tive set of rules are derived by allowing bounded fluctuations in perp space. Furthermore, 
matching rules can be abandoned entirely and quasiperiodicity may arise spontaneously in 
the most probable random tiling [0. 

A fundamental question is whether matching rules are enforced by energetics of real 
materials. Burkov || proposed matching rule enforcement by chemical ordering of Co and 
Cu among certain cites in Al-Co-Cu. However, that model involved an unnatural symmetry 
linking Co sites to Cu sites. Cockayne and Widom || deduced a different, physically realistic, 
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type of Co/Cu ordering based on total energy calculations. In their model, tile edges are 
assigned arrow direction based on their Co/Cu decorations (Fig.^b). The suggested physical 
origin of Co/ Cu chemical ordering rests on the status of Cu as a Noble Metal with completely 
filled d orbitals, unlike normal transition metals such as Co. 

In a tiling model of quasicrystals, the actual atomic interactions in the system Hamilto- 
nian can be replaced with effective interactions between and within tiles HT0| . The resulting 
tile Hamiltonian is a rearrangement of contributions to the actual total energy. In sim- 
ple atomic interaction pictures (pair potentials for example) the relation between the two 
(actual atomic interactions and tile Hamiltonian) is straight forward. It might be difficult 
to find the relations between them for more complicated atomic interactions (many body 
potentials, or full ab-initio energetics, for example) but it is theoretically possible. The tile 
Hamiltonian includes terms which depends only on the number of tiles and other terms 
for different interactions. The tile Hamiltonian greatly simplifies our understanding of the 
relationship between structure and energy, and it is a reasonable way to describe the tiles. 

Space can be tiled in many different ways, even when holding the number of atoms or 
the number of similar tiles fbced. Figure ^ shows three different tiling configurations of 132 
atoms. The first two have the same tiles arranged different ly. The third has the same atoms 
but different tiles. These are called quasicrystal approximants (crystals that are very close 
to quasicrystals in structure and properties). One advantage of approximants is that they 
can be studied using conventional tools developed for ordinary crystals. 



In a previous paper |TT| we studied matching rules in decagonal Al-Co-Cu using a limited 
group of quasicrystal approximants. Some specific details of the tile Hamiltonian couldn't 
be extracted from our limited data set. Here we study more thoroughly the set of rules con- 
trolling these compounds, using different techniques and a much bigger set of approximants. 

We describe our model of decagonal Al-Co-Cu in section y of this paper. Section fT| 
gives our detailed calculations using ab-initio methods. We extract a set of parameters that 
allow an excellent approximation to the total energy. Similar calculations done using pair 



potentials are described in section [IV| for comparison. In section [V|, we talk about various 



other effects that could be considered in a more accurate model. We analyze our findings 
and study their implications for Al-Co-Cu compounds in section |VT|. 



II. DECAGONAL AL-CO-CU MODEL 

Decagonal Al-Co-Cu quasicrystals have been studied by many authors, theoreti- 
cally |8|||[r^] and experimentally fll3"f . Cockayne and Widom || employed mock-ternary 



pair potentials to propose a model based on tiling of space by Hexagon, Boat and Star tiles 
(HBS) decorated deterministically with atoms (Fig. |2|). The tile edge length is 6.38 A. Tile 
vertices are occupied by 11-atom clusters. Each cluster consists of two pentagons of atoms 
stacked on top of each other at |c=2.07 A and rotated by 36°. The pentagon in one layer 
contains only Al atoms. The pentagon in the other layer contains a nrbcture of transition 
metal (TM) atom species and can contain also Al atoms. The mbced Al/TM pentagon 
contains an additional "vertex" Al atom at its center. Ali TM atoms surrounding a vertex 
Al atom belong to tile edges. Decagonal clusters meet along pairs of Co/Cu atoms. It was 
shown in the original model || that TM atoms prefer to alternate in chemical species on tile 
edges. This was confirmed later by ab-initio calculations using the Locally Self-consistent 
Multiple Scattering (LSMS) method [|TT , 14 1 , and is confirmed again in this study. 



We assign arrows to edges based on their TM atom decorations. By our definition, an 
arrow points from the Cu atom towards the Co atom. Tile edges meet in vertices of 72° or 
144° angles. An angle is of type "i" if both edges point in towards their common vertex. 
Types "o" and "m" are out- and mbced-pointing, respectively. 

The HBS tiles are composed of Penrose rhombi with double arrow matching rules sat- 
isfied (by definition) inside the HBS tiles. Some of their properties are summarized in ta- 
ble |. Quasiperiodic tilings can then be constructed from HBS tiles obeying the single-arrow 
matching rules (Fig. 

We choose to define a tile Hamiltonian for the system 
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H = £ + A 72 £ <£? 2 + A 144 Y, nfuEL- (1) 

« /3 /3 

Here is the number of specific tile type a which can be H, B or S; n 72 and T^144 £1X6 the 
numbers of 72° and 144° angles and j3 defines the angle type which can be i, o or m. We fit 
energy parameters in the Hamiltonian (j^) to achieve H ~ E tot for a wide class of structures. 
The tile energies are more important than the other terms, and it turns out that the 72° 
angle interactions are more important than the 144° 's. It is desirable to check how each term 
can alter the system energies. We set A 72 and Ai 44 to 1 or for the purpose of including and 
excluding the 72° and 144° angle interactions. 

An H tile contains 25 atoms, counting those on the perimeter fractionally. Ali but 3 
internal atoms (one Al and two Co) belong to vertex decagonal clusters. The internal Co 
atoms occupy symmetric sites, but the internal Al atom breaks the symmetry by residing 
on one of two geometrically equivalent positions between the two Co atoms (Fig. [I|b, left). 
Its interaction with tile edge TM atoms determines the preferred position. A B tile has 41 
atoms. An internal Al atom breaks the symmetry (Fig. [l]b, center) by residing in one of 
two equivalent sites. An S tile has 57 atoms. Two internal Al atoms break the symmetry 
(Fig. [I|b, right) by occupying any two of five equivalent sites as long they are 144° apart. A 
phason flip can switch a specially arranged star and hexagon into a pair of boats. This is 
shown in Fig. |2]c outlined by a dashed line (compare with Fig. [2]b). 

The present structure differs slightly from the original model ||. In the original model, 
Cu atoms take certain symmetry-breaking positions inside the boat and the star which makes 
it difficult to parameterize their interactions. We choose to replace the Cu atoms in the tile 
interiors with Al atoms. Also in the original model, symmetry-breaking Al atoms inside H, 
B and S were placed in averaged sites between two Co atoms. Their vertical heights lay 
midway between the two main atomic layers. Here we place them in the main atomic layers 
as shown in Fig. |l|b. In terms of the atomic surfaces |J, the atomic surface (AS2) that is 
mainly Al with a thin ring of Cu becomes pure Al, and the Al atomic surface between layers 
(AS3) fills the hole in the pure Al atomic surface (AS2). 
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Many different quasicrystal approximants are exploited here to study different terms in 
the tile Hamiltonian. Ali the approximant tilings we used are listed in table [TT] with some 
of their properties. The smallest tiling is the monoclinic single-hexagon approximant Hi 
(Fig. |3|). It contains one "horizontal" and two "inclined" tile edges. For the decoration 
shown, equation (|IJ) becomes H = E? + \ n (E\ 2 + E° 72 ) + X 1M (2E[ M + 2E° U ). The next 
bigger approximant is a 41-atom single-boat Bi (Fig. ||) for which H = Ef + A7 2 (3.E° 2 ) + 
Ai44(2£'{ 44 + E° iA ) when decorated as shown. Stars alone do not tile the plane, so a single- 
star unit cell approximant is not possible. Two-hexagon approximants can be constructed 
in orthorhombic cells, either by a genuine orthorhombic structure H 2 (Fig. [5]) or by doubling 
the monoclinic cell to create H 2 shown in Fig. |3|. In each case, H = Ef + X 72 (2Ej 2 + 
2E° 2 ) + Ai 44 (4i?^ 44 + 4i?° 44 ). The doubling process gives more freedom in controlling edge 
arrow orientations by performing Co/Cu swaps. Similarly, a two-boat approximant can be 
constructed, either by doubling the single boat tiling Bi or by tiling the two boats as shown 
in B 2 (Fig. |), with H = E? + A 72 (3£ 7 ° 2 + 3E%) + X 1M {3E{ U + E™ 4 + 2E° Ui ) when decorated 
as shown. 

To isolate the tile Hamiltonian parameters Ef and Eff, even larger approximants are 
needed. The three approximants in Fig. each have 132 atoms per unit cell but different tile 
configurations. Two of them (B 2 H' 2 and S1H3 in Fig. |2|b and c respectively) are related to each 
other by a single phason flip. The phason flip turns out to raise the energy, indicating that 
stars are disfavored in these compounds. Angle orientations are investigated by swapping a 
TM pair on tile edges, reversing the directions of the edge arrows. 

Long range interactions and other small terms omitted from the tile Hamiltonian ([]]) can 
be estimated by swapping pairs surrounded by symmetric CoCu bonds on their sides. For 
example, the CoCu bonds on the sides of the horizontal tile edges of H 2 can be specially 
arranged to cancel all angle interactions Eg and leave only other effects. 
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III. AB-INITIO STUDY 



For our calculations we employ ab-initio pseudopotential calculations utilizing the Vienna 
Ab-initio Simulation Package (VASP) program |T^|. We use ultrasoft Vanderbilt type pseu- 
dopotentials [16| as supplied by G. Kresse and J. Hafner [I7||. Our calculations are carried 
out on the Cray T3E and on the newly installed Compaq TCS machine at the Pittsburgh 
Supercomputer Center. 

The k-space mesh size (among other parameters) determines the accuracy of the cal- 
culations. Bigger k-space grids are more accurate but more expensive in calculation time. 
One has to find a balance between the number of atoms in a unit cell and the size of the 
k-space grid in order to fit within the available computer resources. As explained before, 
we use many different approximants for our study each with its own convergence behav- 
ior. The k-space mesh is increased until a consistency of about 0.02 eV is reached in worst 
cases. But within a reasonable use of our allocated computer times we are able to get better 
convergence (0.002-0.01 eV) for most structures. 



Convergence test calculations are summarized in table |ITT| for our Bi tiling and in table [TV 
for our B2H2 tilings. We compared two slightly different structures for each tiling, differing 
in orientation of a single arrow (CoCu pair circled in Fig. |]and Fig. ^a). The tables suggest 
that at the chosen k-space grid used in our calculations, marked by a star in the tables, the 
accuracy is better than 0.02 eV. The smallest grid of 1x1x1 (single k-point in the center of 
the k-space unit cell) for our Bi approximant takes about 8 minutes on the T3E machine 
(450 MHz alpha processors) using 8 processors. The largest grid we use for Bi of 5x5x15 
(188 independent k-point) takes 4.8 hours on 64 processors. The B 2 H 2 structure has 132 
atoms per unit cell. The smallest 1x1x1 grid takes 2 hours on 8 processors. The largest grid 
of 2x2x10 (20 independent k-points) takes about 6 hours on 64 processors. 

For a fixed number of processors, calculation time grows linearly with the number of 
independent k-points. Calculation time decreases linearly with increasing number of pro- 
cessors only up to about 16 processors. Beyond that, the total charging time (number of 
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processorsxelapsed time) increases notably. Large structures and big k-space grids require 
large numbers of processors because they need more memory. 

Note from tables [TTI| and [TV] that the more isotropic the distribution of the k-space grid 
points along k x , k y and k z , the faster the convergence. For our Bi approximant, meshes that 
most isotropically distribute k-space points are 1x1x3 and its multiples, but for finer meshes 
4x4x11 is slightly more isotropic than 4x4x12. For fixed numbers of k x and k y points, the 
total energy converges towards its limiting value as the number of k z points approaches its 
isotropic value. In all our structures, we choose the most isotropic distribution of k-space 
points possible. 

Many rearrangements of edges are performed by swapping TM atom pairs and the dif- 
ferent structure energies are calculated. Most of these are done for the large approximants 
because they give more configurational freedom. In general, pure 72° angle interaction pa- 
rameters (£72) can ke deduced if the bond to be swapped has an equal number of 144° 
connections on each side (since we can arrange for 144° interaction to cancel) and either a 
single 72° connection which can be attached to any side or two 72° connections attached 
to the same side. An example is outlined by a dashed rectangle in Fig. [^a. The bond is 
surrounded by a single 144° on each side. Both point towards the middle bond, making E 444 
on one side, and E^ 4 on the other. One extra 72° on one side makes E° 2 . If the middle pair 
is swapped, the 144° angle interactions will be conserved. The swap merely switches their 
sides, while the E° 2 becomes E^,. The difference in energy between the swapped and the 
basic structures yields the difference E^ — E° 2 = E a f ter — E be f ore = 0.26 eV. Reversing the 
72° outer bond and then making the swap again can give E l 72 — E^\ = 0.35 eV. For pure 144° 
interactions, we use a middle bond surrounded by two 144° angle on one side and a single 
144° angle on the other side. We arrange for any 72° angles to cancel. The pair surrounded 
by a dashed circle in Fig. |2|a is one example that yields E™ u — E° M = 0.06 eV. 

We compute energies for an over-complete set of structures and use a least square fit 
to determine average parameter values. First, we fit with A72 = A144 = in equation ([!]), 
leaving only three adjustable parameters Ef, Ef and Ef . The fit finds the values of our 
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parameters that minimize the root mean square of the difference between the calculated E to t 
and model H values for an ensemble of tilings with different angle orientations. The fitted 
values Ef , Ef and Ef are shown in the second column of table |V[ The fitting is shown in 
Fig. |7|. The graph shows clearly that the three-parameter fit is not adequate because there 
is a significant variations of E tot among different approximants. The main source of this 
variation is angle interactions that contribute to E to t (as calculated by VASP) but not to 
our model when we set A 72 = A 144 = in Eq. (J]]). 

Note that the values of H, B and S are individually meaningless unless compared with 
pure elemental energies. For example, the data quoted does not include arbitrary offsets 
EATOM [15j of each chemical species. However, the difference 2Ef — Ef 1 — Ef = —1.35 eV 



is meaningful because the offsets cancel out (since S1H3 has the same number of atoms of 
each type as B 2 H 2 ). Thus a pair of boats is favored over a hexagon-star pair. 

The five-parameter fit values are shown in the third column of table |V| where we set 
A 72 = 1 to include Ej 2 and set Ai 44 = to exclude -Ef 44 . The 72° angles are all internal to the 
tiles so their energies can't be separated from the tile energies. Only differences in energies 
can be calculated, such as E l 72 — E 72 and E^ — E l 72 , so we set the lowest energy orientation 
E 72 equal to zero without any loss of generality. The fourth column of table [VT] shows the 
eight-parameter fit when we set A72 = A144 = 1. In contrast to the number of 72° angles, the 
number of which is determined entirely by the number of tiles, the number of 144° angles 
depends on the arrangements of tiles. As a result, we may calculate all three energies E^ 
independently. The fitting is shown in Fig. |8|. The remaining deviation from the y=x line 
(standard deviation=0.0013 eV) is due to other effects not included in our model H (Eq. p, 
as well as incomplete convergence or other calculational inaccuracies. 

IV. PAIR POTENTIALS 

The ground state total energy of a system can be expanded in terms of a volume energy 
and potentials describing n-body (n=2,3,4,...) interactions jT3]. The volume energy is the 



9 



dominant contribution to cohesive energy. It depends on the composition and density, but 
not the specific structure. The n-body interactions distinguish between different crystal 
structures at the same composition and density. It is customary to truncate this n-body 
series at the pair potentials (n=2), because they are much easier to calculate and use, and 
because higher order interactions are often weak. Even for transition metals, with their 



localized d-band, the truncation at pair potentials proved to be practical |19|,|20| . Pair 
potentials are functions V a ^{r) of pair separation r and atom types a and (3. 

Many different pair potentials have been used to study this and other quasicrystals. 
Cockayne and Widom U proposed mock-ternary potentials extracted from Al-Co pair po- 
tentials. Noting that, in Al-Co-Cu, Cu substitutes for an equal combination of Al and Co, 
they approximated Cu interactions by the average interactions of Al and Co. In addition, 
the Co-Cu interaction was defined as the average of the Co-Co and Cu-Cu interactions in 
order to obtain ternary potentials from the A1TM binaries. They adopted Al-Co pair po- 



tentials calculated by Phillips et al. ||21|| . Their discovery of alternation of CoCu pairs atoms 
on tile edges, and many other details, are all consistent with our VASP results. 

Later, more rigorous pair potentials derived by Generalized Pseudopotential Theory 
(GPT) were developed for Al-Co- Ni and Al-Co-Cu [T9| , |2"2"| ] . The original GPT pair potentials 
suffered from TM over-binding which is an unphysical attraction between TM atoms at small 
separations. The strongest over-binding appears in Co-Co pair potentials. We modified 
the CoCo and NiNi pair potentials at short distances by adding a repulsive term using 
VASP to get the energy and length scale [pl|. The resulting potentials behave really well 



in simulations [p3[| . The Co-Cu pair potentials were defined as equal to the Ni-Ni pair 
potentials l/ CoCu (r) = V mm( - r \ These, in turn, were close to the average of Co-Co and 
Cu-Cu potentials because Ni resides between Co and Cu in the periodic table. Specifically, 
\/ NiNi (r) w |(y CoCo (r) +y CuCu (r)), with biggest error of 0.002 eV at 3.12 A which is about 
15% error. The other AICoCu pair potentials were found to be well behaved up to large Cu 
composition JEJ. 

We calculated the energies of the approximants using both mock-ternary and modified 
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GPT pair potentials to check how pair potential results compare to VASP. Results of the 
fitting are summarized in table [VT] for all the methods used. Aside from a difference in 
energy scale between the modified GPT and the mock-ternary pair potential calculations, 
they are qualitatively close to each other and to VASP. The order of 144° angle interaction 
is reversed compared to VASP, but these interactions are very weak. 

In table ^ we see that E™ 2 k \[E\ 2 + E° 2 ) and E™ A w \{E\ 4A + E° M ) for all three 
calculation methods (VASP, mGPT, and mock-ternary). To understand this, note that when 
two tile edges meet at a vertex, the TM bonds on them are at three different separations 
from each other. One separation length is between the TM positions close to the vertex, 
r Q is between the far positions and r m is the separation between mbced positions. Take 
the smallest of all, example and consider pair interactions. In bonds with the 

"i" configuration, two Co atoms are distance r j from each other. The energy contribution 
due to this pair is V CoCo (ri). The same positions are occupied by two Cu atoms in the "o" 
configuration with energy contribution l /CuCu (r i ), and by one Co and one Cu atom in the "m" 
configuration with energy V ACuCo (r i ). The contribution to the energy difference E™-\(E\ 2 + 
Ej 2 ) calculated from these pairs at separations n is V CuCo (ri) - \{V CoCo {ri) + V CnCn (ri)). 
Similar identities hold at the separations r Q and r m . In mock-ternary pair potentials, these 
differences of potentials are defined to be zero, sugg esting Ef = + E$) should hold 



exactly. The small deviations from this identity in the fourth column of table |V| are due 
to the fitting procedure. The potential differences are again close to zero for mGPT pair 
potentials as discussed above, and the small deviations from the energy identity in the third 
column of table [V]] are also due to the fitting. In VASP, energies are calculated accurately, 
considering all n-body interactions. An averaging of interactions is not assumed a priori, 
but our calculations confirm that averaging is a good approximation, as shown in the second 
column of table |VT] . 

Another interesting near-degeneracy occurs in a zigzag of 72° angles. Such a zigzag runs 
vertically across Fig. |3|. Consider three consecutive bonds in a zigzag. If both outer bonds 
point in towards, or both point out from, the middle bond, a swap of the middle bond leaves 
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unchanged the total number of "i" , "o " or "m" interactions. If one of the outer bonds point 
in towards the middle bond and the other points out from it, then the swap of the middle 
bond changes the energy by 2E^ — E\ 2 — E% 2 . Because E^ is very close to the average of 
E\ 2 and E° 2 (as previously shown) these configurations are again nearly degenerate. 



V. OTHER EFFECTS 

Chemical ordering of TM atoms on tile edges define edge arrowing in our model. We 
study chemical ordering here using our H 2 approximant (the orthorhombic unit cell in Fig. |3|) 
which contains two horizontal tile edges. When a Cu atom from one horizontal edge is 
swapped with the Co atom on the other, the resulting edges contain pairs of similar species 
(one CoCo and one CuCu). The process raises the energy by 0.68 eV/cell. The same swap 



was studied before with LSMS and gave 0.17 eV/cell. Using the pair potentials, TM 
atoms favor alternation on the tile edges by 0.022 eV/cell for mGPT and 0.079 eV/cell for 
mock-Ternary. Although the magnitude is not certain, the sign consistently favors Co/Cu 
ordering. 

In Al-Co-Ni, CoCo and NiNi pairs are slightly preferred over CoNi pairs. As a result 
Al-Co-Ni has no arrow decorations at low temperatures. Cu and Ni are adjacent in the 
periodic table, but they are notably different in their properties. In an isolated Ni atom, the 
3d shell has six electrons and the 4s shell is filled. The partial filling of the d-band strongly 
influences atomic interactions. The 3d shell in Cu is filled with electrons and the 4s has 
one electron which makes Cu act more like a simple metal. The d-band of Cu is buried 
and doesn't participate strongly in interactions. This is why chemical ordering is strong for 
CoCu pairs but not for CoNi pairs. 

An important issue is the position of the symmetry-breaking Al atom inside a hexagon 
we mentioned in section [TT|. There are two symmetrically related positions between the two 
internal Co atoms, and we force the Al atom to take one of these positions as shown in 
Fig. |I|b (left). If the horizontal edge arrows are parallel to each other, the Al atom prefers 
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to reside in the side closest to the Co atoms by about 0.03 eV. With the off-center Al, we 
define the decomposition of the hexagon into rhombi such that the symmetry-breaking Al is 
placed as in Fig. |1|. The position of the internal Al atom inside a hexagon, together with the 
horizontal tile edge arrows, define a "direction" for the hexagon. We noticed that generally 
hexagons prefer to align parallel to each other in our H 2 structure by about 0.01 eV. These 
effects are very small but are enough to account for some of the discrepancies between E tot 
and H in our calculations. 

The decomposition of the hexagon into rhombi is lost by placing the Al atom exactly at 
the center of the hexagon. However, this position is lower in energy by 0.2 eV as calculated 
by VASP. In the pair potential picture, the central position for the Al atom is preferred 
by 0.11 eV/cell using mGPT and 0.01 eV/cell using mock-ternary potentials. Depending 
on the edge decoration, this Al may relax very slightly from the central positions, but this 
effect is minimal and does not significantly influence the energy. Thus a more realistic model 
in which the Al atoms are centered should be described even more accurately by our tile 
Hamiltonian. 

One more small effect appears in "hidden" 144° angles, where two 72° angles share one 
edge making an extra 144° (Fig. ^ has two hidden 144's). The shared edge orientation affects 
the angle interaction £f 44 of the outer edges. We calculate the difference E" 44 -E^ 4 with the 
shared edge pointing outward and again with it pointing inward. With the shared edge 
outward pointing, the difference E° 44 -E^ 4 = —0.075 eV. An inward-pointing middle edge 
raises the difference by 0.015 eV, so that E° u — E^ = —0.060 eV. This effective three-arrow 
interaction can account for more of the remaining small discrepancies between E tot and H. 

So far we have examined interactions within the quasiperiodic plane. Now consider 
perpendicular interactions. Pairs of TM atoms on tile edges are 1.51 A apart within the 
quasiperiodic plane and 2.07 A apart along the perpendicular, periodic direction. The net 
bond length is 2.56 A. The lines connecting them make a zigzag of alternating TM atoms 
extending along the periodic axes. We turn our attention to atomic order in this direction. 
Approximant B2H2 (Fig. §a) has a horizontal glide plane parallel to the long side of its unit 
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cell that can be exploited for this purpose. We swap one CoCu pair on a horizontal edge 
(call this pair a) and call the structure (A). Another structure (B) is made from B 2 H 2 by 
swapping instead the glide-equivalent image of pair a (call this pair b). These two structures 
have equal energies by symmetry. Further, we build a 264 atom unit cell by stacking two 
132 atom unit cells. It is built once by stacking an A layer over an A layer and another 
time by stacking a B layer over an A layer. In the AA stacking, TM alternation along the 
vertical zigzag is conserved. In AB the zigzag sequence is violated along pair a and along 
pair b. Along each pair the alternation defect includes a CoCo pair and a CuCu pair. The 
AA and AB structure energies are calculated with a k-point mesh of 2x2x5. The difference 
E A b ~ E AA = 0.392 eV per 264-atom cell. 



VI. DISCUSSION 

We discuss here the implications of our fmdings on the structure of decagonal AICoCu. 
The main result is that now energy can be calculated quickly and accurately for these 
compounds by adding the relevant terms in the tile Hamiltonian H (Eq. [1]) using parameters 



obtained in table [VI]. For example, consider the cohesive energy of each tile type. We define 
a tie-line energy E tie _u ne to be the energy per atom of the pure element: fcc Al, fcc Cu and 
spin-polarized hcp Co. The structure energies lie below the tie-line and the difference is the 
cohesive energy per atom, E coh . We calculate £? t(rf [Al] = —4.17 eV/atom, S tot [Cu] = —4.72 
eV/atom and E to t[Co\ = —8.07 eV/atom, all at the experimental lattice constants. The tile 
cohesive energies are i? co h[H] =-7.75 eV , ^^[B] =-12.4 eV and ^^[S] =-16.81 eV (using 
data from our eight-parameter fit). The difference between two boats and a hexagon-star 
pair is 2£' co / l [B] — ^^[H] — -E co h[S] = -0.24 eV. We can add up the cohesive energies of the 
tiles to obtain a quick estimate of the cohesive energy of the quasicrystal. For HBS tilings, 
the "golden" ratio H:B:S=v / 5r:\/5:l can be obtained, for example, by removing double- 



arrow edges from a Penrose tiling |23|.[24|]. For such a tiling the cohesive energy is -0.3035 



eV/atom. Our results show that stars are disfavored, and a tiling with hexagon and boats 
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is lower in energy. The ratio of H:B m HB tilings is l:r and the cohesive energy is -.3045 
eV/atom. 

Most bonds participate in combinations of 144° and 72° angles. The stronger interactions 
determine bond arrowing. When a bond is surrounded by a total of four 144° angles and no 
72° angles, the middle bond is a part of 144° zigzag and its decoration doesn't matter. An 
example of this is circled in Fig. |(| There is only one configuration where a bond orientation 
is determined by 144° interactions. This is the configuration we used to get pure 144° angle 
effects (see sec. These configurations occur occasionally (one is in Fig. ||b), but usually 
bond orientations are determined primarily by 72° interactions. 

Quasicrystals are observed to be stable mainly at high temperatures f25| . This can be 



due to a variety of entropic contributions. Transitions from crystal to quasicrystal phases 
are reported at about T^lOOOK [^fj or about k B T=0.1 eV. At such temperatures the 144° 
angle interactions are irrelevant because they are small compared to energy fluctuations, 
and the structure is determined primarily by its tile types and by the 72° angle interactions. 

Our model expectations are in reasonable agreement with calculated energies, suggesting 
that we have captured the most important energetic effects. The worst deviation is about 
0.1 eV. Out of that we account for 0.03-0.05 eV from the internal Al atom effects on tile 
edges. The rest can be a collection of long range interactions. We do see these long range 
effects in some instances. For example, when calculating pure 72° angle interactions using 
the bond surrounded by a rectangle in Fig. |2|a in two different approximants (B 2 H 2 and 
SiH 3 ). The environments are identical up to about 7 A, but a difference of about 0.02 eV 
in EI^-E^ between the two cases shows up. 

Pair potential calculations show that they are capable of catching qualitatively the dom- 
inant 72° interactions we are investigating with a lot less calculation time. 

In our previous paper |11| , we reported several results related to edge arrowing calculated 



using an all-electron multiple scattering method known as LSMS ||27|| . Approximants H2 and 
H 2 were used, with internal Al atoms centered. The swap energy for chemical ordering agrees 
in sign, but VASP's is four times bigger that LSMS. Other swaps that give 2E r ^ 2 -E l 72 -Ej 2 
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agree in sign, with a similar factor disagreement in magnitudes. 

Further studies might include the effect of TMA1 (as opposed to CoCu) arrows on tile 
edges. The difficulty comes from the fact that such arrowing exist not only on tile edges but 
also inside the tiles. Phason disorder along the periodic axes is important. So far we studied 
only Co/Cu disorder along the periodic axis but not tile flips. The system's behavior under 
relaxation and the preferred relaxed atomic positions are wide areas to explore. Relaxation 
may alter the quantitative values of our tile Hamiltonian parameters. Finally, the biggest 
unresolved question is: what type of structure minimizes the value of our tile Hamiltonian? 

ACKNOWLEDGMENTS 

The authors thank CL. Henley and Y. Wang for useful discussions. IA wishes to thank 
King Abdul Aziz University (Saudi Arabia) for supporting his study and MW acknowledges 
support by the National Science Foundation under grant DMR-01 11198. We thank the 
Pittsburgh Supercomputer Center for computer time used for this study. 



16 



REFERENCES 



[1] P. Gummelt, Geometria Dedicata 62, 1 (1996); P. Gummelt, Material Science and 
Engineering A, 294-296, 250 (2000) 

[2] PJ. Steinhardt and H.-C Jeong, Nature (London) 382, 431 (1996); H.-C Jeong and PJ. 
Steinhardt, Phys. Rev. B 55, 3520 (1999) 

[3] B. Grunbaum and G. C. Shephard, Tilings and patterns (W. H. Freeman, New York , 
1987) 

[4] R. Penrose, Bull Inst. Maths. Its Appl. 10, No. 7/8, 266 (1974) 

[5] M. Gardner, Sci. Am. 236, 110 (1977) 

[6] J.E.S. Socolar, Commun. Math. Phys. 129, 599 (1990) 

[7] CL. Henley, in Quasicrystals: the State of the Art, edited by P.J. Steinhardt and D. P. 
DiVincenzo (World Scientific, Singapore, 1991) 

[8] S. E. Burkov, Phys. Rev. B 47, 12325 (1993) 

[9] E. Cockayne and M. Widom, Phys. Rev. Lett. 81, 598 (1998) 

[10] M. Mihalkovic, W.-J. Zhu, CL. Henley and M. Oxborrow, Phys. Rev. B. 53, 9002, 
(1996) 

[11] M. Widom, I. Al-Lehyani, Y. Wang and E. Cockayne, Matter. Sci. Eng. A295, 8 (2000) 

[12] See for example: A. Yamamoto, K. Kato, T. Shibuya and S. Takeuchi, Phys. Rev. Lett. 
65, 1603 (1990); S. E. Burkov, Phys. Rev. Lett. 67, 614 (1991); M. Torres, G. Pastor, I. 
Jimenez, J. L. Aragon and D. Romeu, Scripta Metallurgica et Materialia 27, 83 (1992); 
G. Trambly de Laissardiere and T. Fujiwara, Phys. Rev. B 50, 9843 (1994); G. Zeger 
and H.-R. Trebin, Phys. Rev. B 54, R720 (1996) 

[13] See for example: L.X. He, Y.K. Wu and K. H. Kuo, J. Mater. Sci. Lett. 7, 1284 (1988); 

17 



A.P. Tsai, A. Inoue and T. Masumoto, JIM Mater. Trans. 30, 300 (1989); H. Chen, 
S. E. Burkov, Y. He, S. J. Poon and G.J. Shiflet, Phys. Rev. Lett. 65, 72 (1990); A.R. 
Kortan, R.S. Becker, F.A. Thiel and H.S. Chen, Phys. Rev. Lett., 64, 200 (1990); W. 
Steurer and K. H. Kuo, Acta Cryst. B46, 703 (1990); K. Hiraga Sci. Rep. RITU A, 36, 
115 (1991); B. Grushko, R. Wittmann and K. Urban, J. Mater. Res. 9, 2899 (1994); K. 
Saitoh, K. Tsuda, M. Tanaka, A.P. Tsai, A. Inoue and T. Masumoto. Phil. Mag. A 73, 
387 (1996); K. Edagawa, K. Suzuki and S. Takeuchi, Phys. Rev. Lett. 85, 1674 (2000) 

[14] For the LSMS method see: Y. Wang, G. M. Stocks, W.A. Shelton, D.M.C. Nicholson, 
Z. Szotek, and W.M. Temmerman, Phys. Rev. Lett, 75 2867 (1995) 

[15] G. Kresse and J. Hafner, Phys. Rev. B 47, RC558 (1993); G. Kresse, Thesis, Technische 
Universitat Wien 1993; G. Kresse and J. Furthmuller, Comput. Mat. Sci. 6, 15-50 
(1996); G. Kresse and J. Furthmuller, Phys. Rev. B 54, 11169 (1996) 

[16] D. Vanderbilt, Phys. Rev. B 41, 7892 (1990) 

[17] G. Kresse and J. Hafner, J. Phys.: Condens. Matter 6, 8245 (1994) 

[18] J. Hafner, From Hamiltonians to Phase Diagrams (Springer-Verlag, Berlin Heidelberg, 
1987) 

[19] I. Al-Lehyani, M. Widom, Y. Wang, N. Moghadam, G. M. Stocks and J.A. Moriarty, 
Phys. Rev. B64, 075109 (2001) 

[20] A. H. MacDonald and Roger Taylor, Can. J. Phys. 62, 796 (1984) 

[21] R. Phillips, J. Zou, A.E. Carlsson, and M. Widom, Phys. Rev. B 49, 9322 (1994) 

[22] J.A. Moriarty and M. Widom, Phys. Rev. B 56, 7905 (1997); M. Widom and J.A. 
Moriarty, Phys Rev. B 58, 8967 (1998); M. Widom, I. Al-Lehyani and J.A. Moriarty, 
Phys Rev. B 62, 3648 (2000) 

[23] M. Mihalkovic, I. Al-Lehyani, E. Cockayne, CL. Henley, N. Moghadam, J.A. Moriarty, 



18 



Y. Wang and M. Widom, Phys. Rev. B 65, 104205 (2002) 
[24] CL. Henley, Phys. Rev. B 34, 797 (1986) 

[25] S. Ritsch, C. Beeli, H.-U. Nissen, T. Godecke, M. Scheffer, and R. Luck, Phil. Mag. 
Lett., 78. 67 (1998) 

[26] Dong C, J. M. Dubois, M. De Boissieu, and C. Janot, Journal of Physics: Condensed 
Matter, 3, 1665 (1991); K. Hiraga, Wei Sun, and F. J. Lincoln, Japanese Journal of 
Applied Physics Letters,30, L302 (1991); K. Hiraga, F.J. Lincoln, and Wei Sun, JIM, 
32, 308 (1991) 

[27] Y. Wang, G. M. Stocks, D. M. C. Nicholson, and W.A. Shelton, Computers, Mathematics 
with Applications, 35,85 (1998) 



19 



TABLES 



TABLE I. Basic tiles in HBS model and their compositions. 



Tile 


Composition 


Penrose Rhombi 


H 


Ali 7 Co 5 Cu 3 


2T+F 


B 


Al 2 9Co 8 Cii4 


T+3F 


S 


AI41C011CU5 


5F 
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TABLE II. The approximant tilings we use for our study, their compositions and the number 
of different decorations (N^) of each of them. The unit cells are either orthorhombic (a, b are given) 
or monoclinic (a, b and 9 are given). Ali of them have c=4.14 A in the z-direction. The number of 
independent k-points is the number on which most of the structures are calculated. To investigate 



the convergence we go higher for a few structures (see tables III and |TV|) . 



Tiling 


Figure 


Composition 


a, b {9) 


Indep. K-points 


N d 


Hi 


i 


AI17C05CU3 


12.14, 7.50 (72°) 


132 


3 


Bi 


1 


AI29C08CU4 


12.13, 12.13 (108°) 


88 


6 


H 2 


i 


Al 34 CoioCu 6 


23.08, 7.56 


66 


5 


H' 2 


1 


Al 3 4Coi Cu 6 


14.27, 12.13 


66 


4 


B 2 


1 


Al 58 Coi 6 Cu 8 


12.13, 30.30 (134.5°) 


55 


3 


B2H2 


|a 


AI92C026CU14 


19.63, 23.08 


20 


27 


B 2 H 2 


0b 


AI92C026CU14 


19.63, 23.08 


20 


6 


S1H3 


|c 


AI92C026CU14 


19.63, 23.08 


20 


12 
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TABLE III. Energies of approximant Bi (Fig. || and one of its single-swap variants). Our 
convergence investigation goes through several k-point grids. Nearly isotropic k-point distributions 
are the 1x1x3 mesh and its multiples. For finer grids, 4x4x11 is more isotropic than 4x4x12. 



K-point grid 


Indep. K-points 


E_Bi 




AE=E^-E Bl 


1x1x1 


1 


-186.83238 


-186.78158 


0.05080 


1x1x2 


1 


-217.02636 


-216.68522 


0.34114 


1x1x3 


2 


-215.95009 


-215.63608 


0.31401 


2x2x2 


4 


-218.18586 


-217.90238 


0.28348 


2x2x4 


8 


-216.75311 


-216.51774 


0.23537 


2x2x6 


12 


-216.74436 


-216.49251 


0.25185 


3x3x3 


14 


-216.59210 


-216.32712 


0.26498 


3x3x9 


41 


-216.74029 


-216.48619 


0.25410 


4x4x4 


32 


-216.77625 


-216.53784 


0.23841 


4x4x11* 


88 


-216.74116 


-216.48698 


0.25418 


4x4x12 


96 


-216.74323 


-216.48616 


0.25707 


5x5x15 


188 


-216.74291 


-216.48615 


0.25676 
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TABLE IV. Energies of approximant B2H2 (Fig. |2|b and one of its single-swap variants.) 



K-point grid 


Indep. K-points 


Eb 2 h 2 


B% H2 


AE=E% w 2ll2 -E B2H2 


1x1x1 


1 


-606.30421 


-606.22773 


0.07648 


1x1x2 


1 


-704.05048 


-703.77620 


0.27428 


1x1x4 


2 


-700.48922 


-700.24548 


0.24374 


2x2x2 


4 


-704.02670 


-703.77286 


0.25384 


1x2x8 


8 


-700.38011 


-700.12094 


0.25917 


2x2x4 


8 


-700.73700 


-700.49654 


0.24046 


1x2x10 


10 


-700.42044 


-700.15409 


0.26635 


2x2x8 


16 


-700.44533 


-700.18672 


0.25861 


2x2x10* 


20 


-700.45159 


-700.19033 


0.26126 
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TABLE V. Fitting our data with different number of parameters. Units are eV. The Standard 
deviation for each set is reported in the last row. Including 72° the 144° angle interactions improved 
the fits as shown. 



Parameters 


A72 — A144 — 


A72 — 1, A144 — 


A72 — A144 — 11 


Ef 


-133.17 


-133.15 


-133.15 


Ef 


-216.42 


-216.76 


-216.77 


Ef 


-298.32 


-300.08 


-300.15 


2Ef-Ef-Ef 


-1.35 


-0.29 


-0.24 


%2 


- 


0.55 


0.55 


%2 




0.23 


0.22 


^72 










F* 
^144 






0.037 


^144 






-0.003 


^144 






-0.034 


Standard deviation 


0.37 


0.0026 


0.0013 
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TABLE VI. The energy costs of significant parameters in our model obtained from 8-parameter 



fit 



Parameters 


Energy (VASP) 


Energy (mGPT) 


Energy (mock-T) 




(eV) 


(eV) 


(eV) 


2Ef-Ef-Ef 


-0.24 


-0.59 


-0.13 




0.55 


0.60 


0.12 


T7 1 T!l 

^72 


0.23 


0.31 


0.05 


^72 











E 144 


0.037 


-0.042 


-0.0079 


H-144 


-0.003 


0.010 


0.0001 


E 144 


-0.034 


0.032 


0.0080 
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FIGURES 




Co • Cu o Al o 



FIG. 1. HBS tiles and their decompositions to Penrose tiles (a) and atomic decorations (b). In 
(b), only TM and symmetry breaking Al atoms are shown. 



26 



27 



FIG. 2. Space can be tiled in many ways using HBS tiles. Ali these approximants (a=E>2H2, 
b=E>2H 2 and c=SiH3) have 132 atoms per unit cell. Structures in (b) and (c) differ by a phason 
flip outlined in (c) with a small dashed line. Bonds surrounded by the dashed square and circle in 
(a) are bonds that can give information about pure angle interactions. 




FIG. 3. Filling space with hexagons. This approximant has one hexagon per monoclinic cell 
(fine dashing, Hi structure in the text). The unit cell has 25 atoms. The cell can be doubled to 
get a 50-atom orthorhombic unit cell (coarse dashing, H2 structure in the text). 
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FIG. 4. Single-boat (Bi) approximant containing 41 atoms per monoclinic cell. The circled 
CoCu pair was used for our convergence study. 
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FIG. 5. Orthorhombic two-hexagon approximant (H2) with a different arrangement of hexagons 
from fig 0. 



30 




FIG. 6. Two-boat approximant (B2). One of the "keel" bonds (circled) in the boats participates 
only in 144° angles and has a highly symmetric Al environment 
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-5.32 -5.3 -5.28 -5.26 -5.24 



H/N (eV/atom) 

FIG. 7. Plots of calculated structure energies vs. our model expectations using only tile energies 
(turning off angle interactions). The spread of energies vertically is due to angle interactions not 
accounted for in the model energy when A 72 = A144 = 0. The diagonal line indicates H = E to t- 
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FIG. 8. Including the angle interactions (setting A72 = A144 = 1) greatly improves the fitting. 
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